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Abstract Suppose Mt is a smooth family of compact connected two di- 
mensional submanifolds of Euclidean space E 3 without boundary varying 
isometrically in their induced Riemannian metrics. Then we show that 
the mean curvature integrals 

/ H t dH 2 

J Mt 

are constant. It is unknown whether there are nontrivial such bendings 
Mt- The estimates also hold for periodic manifolds for which there are 
nontrivial bendings. In addition, our methods work essentially without 
change to show the similar results for submanifolds of H n and S n , to 
wit, if Mt = dX t 

d [ H t dH 2 = -kn - ldV(X t ), 

JMt 

where k = — 1 for H 3 and k — 1 for S 3 . The Euclidean case can be 
viewed as a special case where k — . The rigidity of the mean curvature 
integral can be used to show new rigidity results for isometric embeddings 
and provide new proofs of some well-known results. This, together with 
far-reaching extensions of the results of the present note is done in the 
preprint [6] . Our result should be compared with the well-known formula 
of Herglotz (see [5], also [8] and [2]). 
AMS Classification 53A07, 49Q15 

Keywords Isometric embedding, integral mean curvature, bending, 
varifolds 

1 Introduction 

The underlying idea of this note is the following. Suppose J\f t is a smoothly vary- 
ing family of polyhedral solids having edges {E t (k)} fc , and associated (signed) 
dihedral angles [6 t (k)} k . According to a theorem of Schlafli [7] 

J2\E t (k)\j t 0t(k) = 0. 

k 
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In case edge length is preserved in the family, ie 

||*<*)|-0 
for each time t and each k, then also (product rule) 

±Y,\Et(k)\0t(k)=0. 

k 

Should the cWj's be polyhedral approximations to submanifolds A4t varying 
isometrically, one might regard 

J2\ E t(k)\dt(k) 

k 

as a reasonable approximation to the mean curvature integrals 

/2 



/ 



H t dH z 

'Mt 
and expect 

||*(*)| 

to be small. Hence it is plausible that the mean curvature integrals of the Mt 's 
might be constant. In this note we show that that is indeed the case. 

Examples such as the isometry pictured on page 306 of volume 5 of [8] show 
that the mean curvature integral is not preserved under discrete isometries. 

Two comments are in order. The first is that it is very likely that there are 
no isometric bendings of hyper surf aces. One reason for the existence of the 
current work is to produce a tool for resolving this conjecture (as Herglotz' mean 
curvature variation formula can be used to give a simple proof of Cohn-Vossen's 
theorem on rigidity of convex hypersurfaces). Secondly, the main theorem can 
be viewed as a sort of dual bellows theorem (when the hypersurface in question 
lies in H n or S n ): as the surface is isometrically deformed, the volume of the 
polar dual stays constant. This should be contrasted with the usual bellows 
theorem recently proved by Sabitov, Connelly and Walz [4]. 



2 Terminology and basic facts 

Our object in this section is to set up terminology for a family of manifolds 
varying smoothly through isometries. We consider triangulations of increasing 
fineness varying with the manifolds. To make possible our mean curvature anal- 
ysis we associate integral varifolds with both the manifolds and the polyhedral 
surfaces determined by the triangulations. The mean curvature integral of in- 
terest is identified with (minus two times) the varifold first variation associated 
with the unit normal initial velocity vector field. 
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2.1 Terminology and facts for a static manifold M. 

2.1.1 We suppose that M. C M 3 is a compact connected smooth two di- 
mensional submanifold of M 3 without boundary oriented by a smooth Gauss 
mapping n: M — > S 2 of unit normal vectors. 

2.1.2 H: M. —> M denotes half the sum of principal curvatures in direction n 
at points in M so that Hn is the mean curvature vector field of M . 

2.1.3 We denote by U a suitable neighborhood of Ai in IR 3 in which a smooth 
nearest point retraction mapping p: U — > M is well defined. The smooth signed 
distance function a: U — > M is defined by requiring p = ^(p) + <r(p) n(p(p)) for 
each p. We set 

5 = W: [/ -» M 3 

(so that g\M = n); the vector field g is the initial velocity vector field of the 
deformation 

G t :U^R 3 , G t (p)=p + tg(p) for p G U. 



2.1.4 We denote by 



V = v(M) 



the integral varifold associated with M [1, 3.5]. The first variation distribution 
of V [1, 4.1, 4.2] is representable by integration [1, 4.3] and can be written 

5V = H 2 \-M A(-2fl> 



[1, 4.3.5] so that 

<5F( 5 ) = |w 2 (G t (^)) 



t=o 



g-HndH 2 



M 



■J 



-2 / H dTL 

IM 



here 7Y 2 denotes two dimensional Hausdorff measure in 



2.1.5 By a vertex p in M. we mean any point p in M. By an edge (pq) 
in we mean any (unordered) pair of distinct vertexes p, q in M. which are 
close enough together that there is a unique length minimizing geodesic arc 
\pq\ in M. joining them; in particular (pq) = (qp) . For each edge (pq) we write 
d(pq) = {p, q} and call p a vertex of edge (pq) , etc. We also denote by pq 
the straight line segment in M 3 between p and q, ie the convex hull of p and 
q. By a facet (pqr) in M. we mean any (unordered) triple of distinct vertexes 
p, q, r which are not collinear in M 3 such that {pq) , {qr) , (rp) are edges in 
M. ; in particular, (pqr) = (qpr) = (rpq) , etc. For each facet (pqr) we write 
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d(pqr) = {(pq), (qr), (rp)} and call (pq) an edge of facet {pqr} and also denote 
by pqr the convex hull of p , q , r in M 3 . 

2.1.6 Suppose < r < 1 and < A < 1. By a r, A regular triangulation T 
of Ai of maximum edge length L we mean 

(i) a family T 2 of facets in M , together with 

(ii) the family 71 of all edges of facets in T 2 together with 

(iii) the family 7q of all vertexes of edges in 71 
such that 

(iv) pqr C U for each facet (pqr) in T 2 

(v) M. is partitioned by the family of subsets 



p(pqr ~ (pq UqrU rq)) : (pqr) G 7" 2 }> U |p(pg) ~ {p, <?} : (pg) G T\ 

U I {p} : p G T 



(vi) for facets (pgr) G T 2 we have the uniform nondegeneracy condition: if we 
set u = q — p and v = r — p then 



V — I - — r • U , — r 

U / \u\ 



> t\v\ 



(vii) L = sup { |p - q\ : (pq) G T x } 

(viii) for edges in T\ we have the uniform control on the ratio of lengths: 

inf{|p-g| : (pq) G 71} > XL. 

2.1.7 Fact [3] It is a standard fact about the geometry of smooth subman- 
ifolds that there are < r < 1 and < A < 1 such that for arbitrarily small 
maximum edge lengths L there are r, A regular triangulations of M. of max- 
imum edge length L . We fix such r and A . We hereafter consider only r, A 
regular triangulations T with very small maximum edge length L. Once L is 
small the triangles pqr associated with (pqr) in T 2 are very nearly parallel with 
the tangent plane to M at p. 

2.1.8 Associated with each facet (pqr) in T 2 is the unit normal vector n(pqr) 
to pqr having positive inner product with the normal n(p) to M at p. 
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2.1.9 Associated with each edge (pq) in T\ are exactly two distinct facets 
(pqr) and {pqs} in T 2 . We denote by 

n(pqr) + n(pqs) 



\n(pqr) + n(pgs)| 



n(pq) = 

the average normal vector at pq. 

For each (pq) we further denote by 6(pq) the signed dihedral angle at pq be- 
tween the oriented plane directions of pqr and pqs which is characterized by 
the condition 

2sin ^^MTj n(pq) = V + W 

where 

• V is the unit exterior normal vector to pqr along edge pq, so that, in 
particular, 

V • (p — q) = V • n(pqr) = 0; 

• W is the unit exterior normal vector to pqs along edge pq. 
One checks that 

cos6(pq) = n(pqr) • n(pqs). 
Finally for each (pq) we denote by 

p3 



g{pq) = \p-q\- 1 IgdU 1 

Jpq 



the pq average of g; here TC 1 is one dimensional Hausdorff measure in M 3 . 

2.1.10 Associated with our triangulation T of M. is the polyhedral approxi- 
mation 

N[T] = U{pqf- (pgr) € T 2 } 

and the integral varifold 

V[T}= £ v(mt) = v(AA(T)) 

{pqr)€T 2 

whose first variation distribution is representable by integration 

' ■ (0(jpq)\\ , , 



8V[T]= Yj ^LmA 

(pq)er 1 

[1, 4.3.5] so that 



mm) 



E 

{pq)eT 1 L 



\p-q\ 



2 sin 



n(pg) • g(pg) 
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2.2 Terminology and facts for a flow of manifolds Ait 

2.2.1 As in 2.1.1 we suppose that Ai C M 3 is a compact connected smooth 
two dimensional submanifold of R 3 without boundary oriented by a smooth 
Gauss mapping n: Ai — ► S 2 of unit normal vectors. We suppose additionally 
that ip: (—1, 1) x Ai — > M 3 is a smooth mapping with ip(0, p) = p for each 
p € Ai . For each t we set 

ip[t] = ip(t, ■): Ai -» M 3 and M = </?[t](A4). 

Our principal assumption is that, for each i, the mapping ip\t\: Ai — > A4t is 
an orientation preserving isometric imbedding (of Riemannian manifolds). In 
particular, each A4< C K 3 is a compact connected smooth two dimensional 
submanifold of M 3 without boundary oriented by a smooth Gauss mapping 
n t : Ait — > S 2 of unit normal vectors. 

2.2.2 As in 2.1.2, for each i, we denote by H t n t the mean curvature vector 
field of M t ■ 

2.2.3 As in 2.1.3, for each t we denote by U t a suitable neighborhood of Ait 
in M 3 in which a smooth nearest point retraction mapping p t : U t —> Ai t is well 
defined together with smooth signed distance function cr t : {7 t — > M; also we set 
<7[i] = Vcr t : CT* — > M 3 as an initial velocity vector field. 

2.2.4 By a convenient abuse of notation we assume that we can define a smooth 
map 

<p: (-1, 1) x U — ► R 3 , 
p) = p (p) + o"oO)n (p(p)) = po(p)) + <7o(p)n t (p (p)) 

for each t and p. With <^[t] = <p(t, •) we have <p[0] = lu and, additionally, 
cr (p) = a t ((p[t](p)) . We further assume that 

U t = <p[t] U Q 

for each t. 

2.2.5 Fact If we replace our initial <p[i]: — > M 3 's by v 3 !/^^] f° r large enough 
\x (equivalently, restrict times t to —l/fi<t< and decrease the size of Uo 
then the extended tp[t\: Uo — > M 3 's will exist. Such restrictions do not matter 
in the proof of our main assertion, since it is local in time and requires only 
small neighborhoods of the Ait 's. 

2.1.6 As in 2.1.4, for each t we denote by 

V t = v(M t ) 
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the integral varifold associated with M. t . 



2.2.7 We fix < r < 1/2 and < A < 1/2 as in 2.1.7 and fix 2r, 2A 
regular triangulations T(l), 7~(2), 7~(3), ... of M. having maximum edge 
lengths L(l), L(2), L(3) ... respectively with lim^oo L(j) = 0. For each j, 
the vertexes of T(j) are denoted T (j), the edges are denoted 7[(j), and the 
facets are denoted T 2 (j). For all large j and each t we have triangulations 
T(l, t), T(2, t), T(3, i), ... of Mt as follows. With notation similar to that 
above we specify, for each j and t, 



Tot?, = \ <p[t](p) ■■ P G T (j)|, 71 (j, t) = |<¥»[t](p) *>[<]((?)) : (pg> G Mj) 
T 2 (j, t) = (<¥»[t](p)^](?)¥>[t](r)> : (pgr) € T 2 (j) 



2.2.8 Fact If we replace y>[i] by (p[fii\ for large enough u. (equivalently, re- 
strict times t to -1/// < i < 1/fj,) then 7~(1, t), T(2, i), T(3, i), ... will 
a sequence of r, A regular triangulations of Ai with maximum edge lengths 
L(j, t) converging to uniformly in time t as j — > oo. Such restrictions do not 
matter in the proof of our main assertion, since it is local in time. We assume 
this has been done, if necessary, and that each of the triangulations T(j, t) is 
t, A regular with maximum edge lengths L(j, t) converging to as indicated. 



2.2.9 As in 2.1.8 we associate with each j , t, and (pqr) G T 2 (j) a unit normal 
vector n[t, j](pqr) to (f[t](p) f[t](q) <p[t](r) . As in 2.1.9 we assoc iate with each 
j, t, and (pq) £ an average normal vector n[i, j](pq) at <^[t](g) 

and a signed dihedral angle 0[t, j](pq) at y>[i](p) <p[t](g) and the ip[t](p) ip[t](q) 
average g[t, j](pq) of g[t}. 



2.2.10 As in 2.1.10 we associate with each triangulation T(J, t) of M t a 
polyhedral approximation M[T(j, t)] and an integral varifold 



V[T(j, t)} = v{Af[T(j, t)}) = E v ( ^M(P) 

( P gr>er 1 (j) 



with first variation distribution 

<5V[T(j,i)] = £ 

(pq)er 1 (j) 



<p[t])p)<p[t](q) 



A 



2 sin 



n[t, j](pq). 
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SO that 

6V[T(j, t)]( 9 [t]) 

= E [M*](p)-rt*i(«)l 



2 sin 



n[t, j](pq) -g[t, j](pq) 



2.2.11 The quantity we wish to show is constant in time is 

T 



H t dH 2 



Mt 



SV t {g[t}). 



Since, for each time t, 
we know, for each t 



V t = lim V[T(j, t)] (as varifolds) 

j— >oo 



5V t (g[t})= lim 5V[T(j, t)](g[t}). 

We are thus led to seek to estimate 

d 



dt 



5V[T(j, t)](g[t]) 



using the formula in 2.2.10. A key equality it provided by Schlafh's theorem 
mentioned above which, in the present terminology, asserts for each j and t, 

E \\mp)-mq)\ ] 1 

(pq)er 1 (j) 



dt 



0[t,j](pq) 



0. 



2.2.12 Fact Since, for each (ppq) in T 2 (j), d{pqr) consists of exactly three 
edges, and, for each (pq) in ^(j), there are exactly two distinct facets (pqr) 
in T 2 (j) for which (pq) G d(pqr) we infer that, for each j , 

3 

card[Ti(j)] = -card[T 2 (j)] . 

We then use the r, A regularity of the the T(j) 's to check that that, for each 
time t and each (ppq) in T 2 (j) the following four numbers have bounded ratios 
(independent of j, t, and (ppq)) with each other 



H 2 ( <p[t](p) <p[t](q) <p[t](r) ) , \<p[t](p) ~ <p[t](q)\ 
Since 



L(j, t) 2 , L(j) 2 . 

lim n 2 (M[j, t]) = n 2 (M t ) = n 2 (M), 



we infer 



E £(i) 2 <°o, lim E L (i) 3 = - 



sup 

j (pq)eTi(j) 



(pg)er 1 (j) 
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3 Modifications of the flow 



3.1 Justification for computing with modified flows 



As indicated in 2.2, we wish to estimate the time derivatives of 
SV[T(j, t)](g[t]) 



E 

( P q)eT 1 (j) 



\<p[t]<P) - <p[t](Q)\ 



2 sin 



n[i, j](PQ) ■ 9% j](pq) 



In each of the (pq) summands, each of the three factors 



\v[t](p) - <p[t](Q)\ 



2 sin 



n[t, j](pq) -g[t, j](pq) 



is an intrinsic geometric quantity (at each time) whose value does not change 
under isometries of the ambient M 3 . With (pqr) and (pqs) denoting the two 
facets sharing edge {pq) , we infer that each of the factors depends at most on 
the relative positions of (p{t](p) , (p[t](q), ip[t] (r), <p[t](s) and <p[t].M. Suppose 
ip: (—1, l)xR 3 ^ M 3 is continuously differentiable, and for each t, the function 
ip[t] = ip(t, •): M. 3 — > M 3 is an isometry. Suppose further, we set 

<p*(t,p) = 4>{t,<p(t,p)), <p*[t]=<p*(t, •) 

for each t and p so that ip*[t] = ip[t] o ip[t] . If we replace M by M* = ip[0]M 
and (f by ip* then we could follow the procedures of 2.1 and 2.2 to construct 
triangulations and polyhedral approximations T*[j, t] and varifolds V* , etc. 
with 

8V[T(j, t)]{g[t])=SV*[T*(j, t)]{g*[t}). 

Not only do we have equality in the sum, but, for each (pq) the corresponding 
summands are identical numerically. Hence, in evaluating 5V[T(j, t)](g[tf) we 
are free to (and will) use a different tp and ip* for each summand. 



3.2 Conventions for derivatives 



Suppose W is an open subset of M. M and / = (J 1 , f 2 , 
K times continuously differentiable. We denote by 



-.,/"): W 



is 



K 



\\\D 

the supremum of the partial derivatives 

d k f K 



-ip) 



dx i( i)dx i{ 2) ...dxi( K ) 

corresponding to all points p G W , all (i(l), i(2), . . . , i(K)} C {l, . . . , M} 
and k = 1, . . . , N , all choices of orthonormal coordinates {x\, . . . , im) for 
and all choices of orthonormal coordinates (yi, . . . , y^) for 



!>A ; 
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3.3 Conventions for inequalities 

In making various estimates we will use use the largest edge length of the 
jth triangulation, typically called L, and a general purpose constant C . The 
constant C will have different values in different contexts (even in the same 
formula). What is implied is that, with M and ip fixed, the constants C can 
be chosen independent of the level of triangulation (once it is fine enough) and 
independent of time t and independent of the various modifications of our flow 
which are used in obtaining our estimates. As a representative example of our 
terminology, the expression 

A = B ± CL 2 

means 

-CL 2 <A-B< CL 2 . 

3.4 Fixing a vertex at the origin 

Suppose p is a vertex in M. and 

¥>*(— 1, 1) x Uq —> M 3 , <p*(t, q) = <p(t, q) - <p(t, p) for each q. 
Then tp*(t, p) = (0, 0, 0) for each t. One checks, for K = 0, 1, 2, 3 that 

\\\D K (p* HI < 2|||D*' ¥ >|||, HI III = \\\D K ^[t] HI 

for each t. 



3.5 Mapping a frame to the basis vectors 



Suppose (0, 0, 0) G M and that ei and e 2 are tangent to M at (0, 0, 0). 
Suppose also ip(t, 0, 0, 0) = (0, 0, 0) for each t. Then the mapping ip* given 
by setting 

0,0,0) g(t, 0,0,0) ^(t,0,0,0)\ 



v *[t] 



§£(t, 0,0,0) g(t, 0,0,0) |g(t, 0,0,0) 
fet, 0,0,0) 0,0,0) |£(t,0,0,0)/ 



satisfies 
with 



ip* [t] (0, 0, 0) = (0, 0, 0) , Dip* [t] (0,0,0) = 1 R8 



\D K ip* 



\D K v 



for each K = 1, 2, 3 and each t, and 



9t 



<3(iiiDViii-iii^viii + iiPV[*]in 2 )- 
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3.6 Theorem There is C < oo such that the following is true for all suffi- 
ciently small 5 > 0. Suppose 70: [0, 6] — ► M is an arc length parametrization 
of a length minimizing geodesic in M and set 

7(5, t) = <p[i](7o(s)) for each s and t 

so that s — > 7(5, t) is an arc length parametrization of a geodesic in M. t . We 
also set 

r(s, t) = (7(0, t) — 7(5, t)\ for each s and t 
and, for (fixed) < R < 5, consider 

r(R, t) = | 7 (0, t) - 7(12, t)\ for each t. 

Then 



and 



Proof We will show 



^r{R, t) = ±CR 2 



limR-^rfR, t) = 0. 

Rio dt y ' 



= ±CR 2 . 

t=o 



Step 1 Replacing <p(t, p) by ip*(t, p) = tp(t, p) — tp(t, 70 (0)) as in 3.4 if nec- 
essary we assume without loss of generality that 7(0, t) = (0, 0, 0) for each 
t. 

Step 2 Rotating coordinates if necessary we assume without loss of generality 
that ei and e 2 are tangent to A4o at (0, 0, 0) and that 7q(0) = ei 

Step 3 Rotating coordinates as time changes as in 3.5 if necessary we assume 
without loss of generality that Dip[t](0, 0, 0) = 1r3 for each t. 

Step 4 We define 

X(s,t) = j(s,t)-e 1 , Y(s,t) = j(s,t) -e 2 , Z(s, t) = j(s, t) ■ e 3 
so that 

7 (s, t) = (X(s, t), Y(s, t), Z(s, t)) 
and estimate for each s and t: 
(a) X(0, t) = Y(0, t) = Z(0, t) = (by step 1) 
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(b) X t (0, 0) = Y t (0, 0) = Z t (0, 0) = 

(c) X s (s, t) 2 + Y s (s, t) 2 + Z s (s, t) 2 = 1 

(d) X s (s, t) = ±1, y(a, t) = ±1, Z s (a, t) = ±1 

(e) 1/2 < r(s, i)/|s| < 1 (since S is small) 

(f) X(s, 0) = ±Ca, y(s, 0) = ±Cs, Z(s, 0) = ±Cs 

(g) X s (0, t) = X a (0, 0), y (0, t) = Y s (0, 0), Z s (0, i) = Z s (0, 0) (by step 3) 

(h) X st (0, 0) = Y st (0, 0) = Z st (0, 0) = 

(i) X st (s, 0) =X st (0, 0)+ / X sst (r?, 0)^ = 0±ssup|X sst | = ±Ca, 

■/ o 

nt(a, 0) = ±Cs, Z st {s, 0) = ±Cs 

(j) X t (a, 0) = X t (0, 0) + / X st (r?, 0) dr/ = ± Cs 2 , 

Jo 

Y t (s, 0) = ±Cs 2 , Z t (s, 0) = ±Cs 2 

(k) r 2 = x 2 + y 2 + z 2 

{I) rr s = XX S + YY S + ZZ S , r s = ^ (XX S + YY S + ZZ S ) 

(m) rr t = XX t + YY t + ZZ t , r t = ± (XX t + YY t + ZZ t ) 

(n) r s r t + rr st = X s X t + II s( + Y s Y t + yy st + Z s Z t + ZZ st 
(o) evaluating (n) at t = 0, r > we see 

-J^((±Cs)(±l))((±C S )(±C S 2 )) +r(s, 0)r st (a, 0) 

=(±l)(±Cs 2 ) + (±Cs)(±Cs) 

(p) r st (s, 0)=±Cs 

/•it /.fl 

(q) r t (R, 0) =r t (0, 0) + / r st (s, 0) ds = + ±Csds = ±CR 2 . 

Jo Jo 



Geometry & Topology Monographs, Volume 1 (1998) 



The mean curvature integral is invariant under bending 



13 



3.7 Corollary Suppose triangulation T(j) has maximum edge length L = 
L(j) and (pq) is an edge in Ti(j). Then, for each t, 



<p[t](p)-<p[t](q) 



±CL 



and — 

dt 



¥>[t](p) -¥>[*](?) 



±CL 2 



3.8 Stabilizing the facets of an edge 

Suppose T(j) is a triangulation with maximum edge length L = L(j) and that 
(ABC), (ACD) are facets in T 2 (j) as illustrated 

D = (e,f,0) 

/ \ 
(0,0,0)= A a — ► C = (d,0,0). 

X / 
B = (a, b, c) 

Interchanging B and D if necessary we assume without loss of generality the 
the average normal n[0, AC] to M.q at A has positive inner product with 
(C-A)x(D-A). 

1) Fixing A at the origin Modifying cp if necessary as in 3.4 if necessary 
we can assume without loss of generality that <^[i](.A) = (0, 0, 0) for each t. As 
indicated there, various derivative bounds are increased by, at most, a controlled 
amount. 

2) Convenient rotations We set u(t) = <p[t](C), v(t) = ip[t](D) and use 
the Gramm-Schmidt orthonormalization process to construct 

u(t) = ^- v(t) = v(t) ~ v(t) • m m Wit) = U(t) x V(t) 

U[) koi' V{) \v(t)- v (t)-u(t)u(t)y W[) u w* v w- 

One uses the mean value theorem in checking 

III^WIII IIPVIH j , etc 

for each K = 0, 1, 2. We denote by Q(t) the orthogonal matrices having 
columns equal to U(t), V(t), W(t) respectively (which is the inverse matrix to 
its transpose). Replacing ip t by Q(t) o (p t if necessary, we assume without loss 
of generality that there are functions a(t) , b(t) , c(t) , d(t) , e(t) , f(t) , such that 

<p[t](A) = (0, 0, 0), <p[t](B) = (o(t), b(t), c(t)), 

<p[t]{C) = (d(t), 0, 0), <p[t](D) = (e(t), f(t), 0). 
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We assume without loss of generality the existence of functions -F[i](x, y) de- 
fined for (x, y) near (0, 0) such that, near (0, 0, 0) our manifold M. t is the 
graph of F[t] . In particular, 

c(t)=F[t](a{t), b(t)). 

We assert that if \p\ < CL, then 

\F[t](p)\ < CL 2 , \VF[t\(p)\<CL. (3.8.1) 

To see this, first we note that F[t](A) = F[t](C) = F[t](D) = 0. Next we 
invoke Rolle's theorem to conclude the existence of c\ on segment AD and c 2 
on segment CD such 

° 1 , DF[t]( Cl )) = = (£f^, DF[t](c 2 ) 



{D-Ay 7 \\D-C 

Since \p\ < CL we infer 

° A -, DF[t](p)) = ±CL, (^r 1 ^, DF[t](p) ) = ±CL. 



\D-Ay LJV ^7 ' \|D-C| ; 

In view of 2.1.6(vi)(vii)(viii) and 2.2.7 we infer that ei and e 2 are bounded 
linear combinations of (D — A)/\D — A\ and (D — C)/\D — C\ from which we 
conclude that \VF[t](p)\ < CL. This in turn implies that \F[t](p)\ < CL 2 as 
asserted. 

Since 

J^M(o, o) = o 

^F[t\(p) = ±CL (3.8.2) 
|(^](A).e 3 )=0 

c'it) = ^F[t](a(t), b(t)) = | (<p[t](B) • e 3 ) = ±CL. (3.8.3) 

3.9 Proposition Let L, A, B, C, D, a, b, c, d, e, f be as in 3.8. Then 

(1) a'(t) = ±CL 2 

(2) b'{t) = ±CL 2 

(3) c'{t) = ±CL 

(4) d'(t) = ±CL 2 

(5) e'{t) = ±CL 2 

(6) f(t) = ±CL 2 . 



we infer 



and since 



we infer 
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Proof According to 3.7, if r(t) denotes the distance between the endpoints of 
an edge of arc length L at time t, then 

r'(t) = ±CL 2 . 

(i) We invoke 3.7 directly to infer (4) above. 

(ii) We apply 3.7 to the distance between (0, 0, 0) and (e, /, 0) to infer 

i (e 2 + f 2 Y = {ee ' + ff '} = ±CL\ ee> + //' = ±CL 3 . 
dt (e 2 + / 2 ) 5 

(iii) We apply 3.7 to the distance between (d, 0, 0) and (e, /, 0) to infer 

d _ df + n i je-w-*) + f r) = ±CL! _ 
dt {{e-dy + py 

(e-d)(e' -d') + ff = ±CL 3 . 
We subtract the first inequality from the second to infer 

ed' -de' + dd' = ±CL 3 , de' ±CL 3 , e' = ±CL 2 . 
Assertions (5) and (6) follow readily. 

(iv) We apply 3.7 to the distance between (0, 0, 0) and (a, b, c) to infer 

d , n o ox i faa' + bb' + cc') n . , 

— (a 2 + b 2 + c 2 ) 2 = ^ 1 = ±CL 2 , aa' + bb' + cc' = ±CL 3 . 

dt (a 2 + 6 2 +c 2 ) 5 

(v) We apply 3.7 to the distance between (d, 0, 0) and (a, b, c) to infer 

± ((a -d) 2 + b 2 + c 2 Y= (^-d)W-d') + bb' + cc') = ±cl2 ^ 
dt {{a-d) 2 +b 2 + c 2 Y 

(a - d){a' - df) + bb' + cc' = ±CL 3 . 
We subtract the first inequality form the second to infer 

ad' -da' + dd' = ±CL 3 , da'±CL 3 , a' = ±CL 2 , 
which gives assertion (1). 

(vi) We estimate from 3.8 that 

c = F[t](a, b) = ±CL 2 , c' = ^F[t](a, b) + VF[t]{a, b) ■ {a', b') = ±CL, 

which gives (3) above. We have also cc' = ±CL 3 . We recall (iv) above and 
estimate 

aa' + bb' + cc' = ±CL 3 , bb' = ±CL 3 , b' = ±CL 2 , 
which is (2) above. □ 
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3.10 Proposition Suppose T(J) is a triangulation with maximum edge length 
L = L(j) and (pq) is an edge in Ti(j) . Abbreviate 0(t) = 6[t, j](pq) . Then, 
for each t , 



(1) 
(2) 
(3) 
(4) 



9(t) = ±CL 



±CL 



V{t) = ±c 



d_ 
It 



2 sin 



6{t) 



±C 



(5) 



d_ 

dt 



2 sin 



6(t) 



= ±CL 2 



Proof Making the modifications of 3.8 if necessary, we assume without loss 
of generality (in the terminology there) that (p[t](p) = A = (0, 0, 0), y[i](<z) = 
C = (d(t), 0, 0), and that there are (pqB*), (pqD*) 6 T 2 {j)o with </?[£](£?*) = 
B = (a(t), b(t), c(t)), ip[t](D.) = D = (e(t), f(t), 0). 



The unit normal to ACD is (0, 0, 1) while the unit normal to ABC is 

(0, -c, b) 



+ C 2 )5 



so that 



cos# 



(6 2 + c 2 )" 



sin# = ± (1 - cos 2 9) 2 = ± (1 - 



b 2 



±- 



& 2 + c 2 y (6 2 + c 2 

in view of 3.8. Assertions (1) and (2) follow. We compute further 

(b 2 +c 2 )^c' -c bb ' +cc ' 



- = ±CL 



(smO)' = cos 6 6' = ±- 



(6 2 + C 2 )2 



b 2 +C 2 



±c 



in view of 3.9(1) (2) (3) and 3.8. Assertion (3) and (4) follow. Assertion (5) 
follows from differentiation and assertions (1) and (3). □ 
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3.11 Proposition Suppose T(J) is a triangulation with maximum edge length 
L = L(j) and (pq) is an edge in 7i(j). Then 



(1) n[t, j](pq) = {0, ±CL, 1±CL 4 ) 

(2) {d/dt)(n[t, j]{pq)) = (0, ±C, ±CL) + (±CL, ±CL, ±CL) 

(3) g[t, j]( P q) = ( ±CL, ±CL, 1 ± CL 2 ) 

(4) (d/dt)g[t, j](pq) = ( ± C, ±C, 0) + ( ± CL, ±CL, ±CL) 

(5) n[t,j](pq)-g[t,j}(pq) = l±CL 2 



(7) l-n[t,j](pq)-g[t,j](pq)=±CL 2 . 

Proof We let A, B, C, D, F[t], b(t), c(t), d(t) be as in 3.8. We abbreviate 
n = n[t, j](pq) and estimate 



The first assertion follows from 3.8.1. We differentiate to conclude n' = 
±CL(0, -d, V ± C(W + cd)/L - (L/L)(W + cd ± b'L + ±C(b/L){bb' + cd) 

= (0, ±C, ±CL) + (±CL, ±CL, ±CL) 
in view of 3.9(2) (3). This is assertion (2). 
We abbreviate g = g[t, j](pq) and estimate 



(6) (d/dt) n[t, j](pq) ■ g[t, j](pq) = ±CL 



n = 



(0, 0, l) + (0, -c, 6)/(6 2 + c 2 )i 
|(0, 0, l) + (0, -c, b)/(b 2 + c 2 ) l 2 
(0, -c, b+{b 2 + c 2 )^) 



2h (b 2 + c 2 + b{b 2 + c 2 )^y ' 



9 = 



d(t)J |( - F[t] x , -F[t] y , lj| 

1 f d{t) ( - F[t] X , -F[t]y, l) 



1 /■"(*) (-F[t] a ,-F[t] w , 1) 
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The third assertion follows from 3.8.1. We differentiate to estimate that dg/dt 
equals 

-d!_ r d ^ ( - F[t] x , -F[t] y , 1) d' ( - F[t] x , -F[t] y , 1) 

* ^ (1 + + F[t]*) * + d (1 + + F[t]2) * 
1 /" d ±CL(-F[t] te , -F[t] ty ,0) 
+ di l + F[t]l+F[t\l 

1 f d ( -F[t] x , -F[t] y , i)(±C/L)(F[t] x F[t] tx + F[t] v F[t] tv ) 
dj l + F[t]l + F[t]l 

L( ± C, ±C, ±C) + L( ± C, ±C, ±C) + ( ± C, ±C, 0) + L{ ± C, ±C, ±C) 

which gives assertion (4). Assertion (5) follows from assertions (1) and (3). 
Assertion (6) follows from assertions (1), (2), (3), (4) and integration by parts. 
Assertion (7) follows from assertions (1) and (3). □ 



4 Constancy of the mean curvature integral 



4.1 The derivative estimates 



Suppose triangulation T(j) has maximum edge length L = L(j). We recall 
from 2.2.10 that 



5V[T(j, t)]{g[t}) 

= E [m*](p)-vM(?)I 

and we estimate, for each t that 



2 sin 



0[t, j](pq) 



n[t, j](pq) ■ g[t, j](pq) 



d 



j t [6V[T(j) t ](g[t}) 

= E [m*](p)-vM(9)I 

{ P q)er 1 {j) 

+ E fkW(p)-vW(9)l 

(pq)er 1 (j) 

+ E fkW(p)-vW(9)l 

{ P q)er 1 {j) 



2 8in i fiM) 



2 sin 



2 sin 



n [*, j](pq) -9[t, j](pq) 

n[t, j](pq)-g[t, j](pq) 
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We assert that 



d_ 

lit 



(sV[T(j,t)](g[t])) = ±CL3 = E ±^(J) 3 - 

v J { P q)eT 1 (j) { P q)er 1 (j) 

To see this we will estimate each of the three summands above. 

First summand We use 3.7, 3.10(2), 3.11(5) to estimate for each pq, 

[t, j](pq) \ 



Mt](p) -¥>[*](?)! 



2 sin 



2 J 
{CL 2 ){CL){l±CL 2 ). 



j\im) ■ 9[t, jKpq) 



Second summand We use 3.10(5), 3.11(7) to estimate for each pq, 



\<p[t](p) - <p[t]{Q)\ 

<p[t}(p) - <p[t](q)\ 

+ \<p[t](p)-<p[t](q)\ 

+ \<p[t](p)-<p[t](q)\ 

\<p[t](p) - <p[t]{Q)\ 



n[t, j](pq)-g[t, j](pq) 



0% j](pq) 



2sin(«)^[t ljlW 



0[t,j](pq) 



n[t, j](pq) -g[t, j](pq) 
±(CL)(CL 2 )±(CL)(C)(CL 2 ). 



Third summand We use 3.10(2) and 3.11(6) to estimate 



\v\t\ip) - v[t\(q)\ 



j](pq) -g[t, j](pq) 



= (CL)(CL)(CL). 
According to Schlafli's formula [7], 



E 

(pq)eT 1 (j) L 



\f[t]{p) - <p[t]{q)\ 



0[t, j](pq) 



= 0. 



Our assertion follows. 



Qeometry & Topology Monographs, Volume 1 (1998) 



20 



Frederic J Almgren Jr and Igor Rivin 



4.2 Main Theorem 

(1) For each fixed time t, 

lim 5V[T(j, t)](g[t})=5V t (g[t}). 

(2) For each fixed j , 5V[T(j) t ] (g[t]) is a differentiate function of t and 

i lim|(^[r( A ]( 9[t ]))=0 

uniformly in t . 

(3) For each t 

[ H t dH 2 = [ HdH 2 . 
This is the main result of this note. 

Proof To prove the first assertion, we check that 

(fH) t V[T(j, t)\ = V t 

for each t and all large j . Indeed, the r regularity of our triangulations implies 
that the normal directions of the M[T(j) t ] are very nearly equal to the normal 
directions of nearby points on M t and that the restriction of Dp t to the tangent 
planes of the J\f[T(j) t ] is very nearly an orthogonal injection. The first assertion 
follows with use of the first variation formula given in [14.1, 4.2]. Assertion (2) 
follows from 4.1 since 

(pq)eT 1 (j) 

is dominated by the area of A4 (see 2.2.12) and lim^oo L{j) = 0. Assertion 
(3) follows from assertions (1) and (2) and our observation in 2.1.4. □ 

Acknowledgements Fred Almgren tragically passed away shortly after this 
note was written. Since then, the main result for smooth surfaces has been 
reproved in an easier way and generalized to the setting of Einstein manifolds 
by J-M Schlenker together with the second author of the current paper [6]. 
Nonetheless, it seems clear that the methods used here can be used to extend 
these results in other directions. 
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